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We construct a two-dimensional microscopic model of interacting quantum dimers that displays 
an infinite number of periodic striped phases in its T = phase diagram. The phases form an 
incomplete devil's staircase and the period becomes arbitrarily large as the staircase is traversed. 
The Hamiltonian has purely short-range interactions, does not break any symmetries of the under- 
lying square lattice, and is generic in that it does not involve the fine-tuning of a large number 
of parameters. Our model, a quantum mechanical analog of the Pokrovsky-Talapov model of fluc- 
tuating domain walls in two dimensional classical statistical mechanics, provides a mechanism by 
which striped phases with periods large compared to the lattice spacing can, in principle, form in 
frustrated quantum magnetic systems with only short-ranged interactions and no explicitly broken 
symmetries. 



I. INTRODUCTION 

The past two decades have seen the discovery of a 
number of strongly correlated materials with unconven- 
tional physical properties. Due to the competing effects 
of essentially electronic processes and interactions these 
doped Mott insulators typically exhibit complex phase 
diagrams which include antiferromagnetic phases, gen- 
erally incommensurate charge-ordered phases, and high 
temperature superconducting phases. When conduct- 
ing, these systems do not have well defined electron-like 
quasiparticles and their metallic states thus cannot be ex- 
plained by the conventional theory of metals, the Landau 
theory of the Fermi liquid, and the associated supercon- 
ducting states cannot be described in terms of the BCS 
mechanism for superconductivity. 

The startling properties of these materials have led 
to a number of proposals of non-trivial ground states of 
strongly correlated systems which share the common fea- 
ture that they cannot be adiabatically obtained from the 
physics of non- interacting electrons. A class of proposed 
ground states are the resonating valence bond (RVB) spin 
liquid phases, quantum liquid ground states in which 
there is no long range spin order of any kind, and the 
related valence bond crystal phases, of frustrated quan- 
tum antiferromagnetsi^ and their descendants^^ On 
the other hand, the presence of competing spatially- 
inhomogeneous charge-ordered phases in close proximity 
to both antiferromagnetism and high T c superconductiv- 
ity, and the existence of incommensurate low energy fluc- 
tuations in the latter phase, strongly suggest that these 
phases may have a common origin. It has long been sug- 
gested that some form of frustration of the charge degrees 
of freedom may be at work in these systems.—^ The ex- 
planation of both the existence of a large pairing scale in 
the superconducting phase and their close proximity to 
inhomogeneous charge-ordered phases is one of the cen- 
tral conceptual challenges in the physics of these doped 
Mott insulators*^ 

The most studied class of these strongly correlated ma- 



terials are the cuprate high temperature superconductors 
(for a recent review on their behavior and open ques- 
tions see Ref. [nj.) Unconventional behaviors have also 
been seen in other strongly correlated complex oxides^ 
More recently, strong evidence for non-magnetic phases 
has been discovered in new frustrated quantum mag- 
netic materials, including the quasi-2D triangular anti- 
ferromagnetic insulators such as^ CS2CUCI4, the quasi- 
2D triangular organic compounds such asi^ k-(BEDT- 
TTF)2Cu2(CN)3, and the 3D pyrochlore antiferromag- 
nets such as the spin-ice compoun d 15 ' 16 Dy2Ti2C>7 (al- 
though quantum effects do not appear to be prominent 
in spin-ice systems.) 

It is thus of interest to develop a theoretical framework 
to describe quantum frustrated systems in the regime of 
strong correlation, and to understand their role in the 
mechanism for inhomogeneous phases in strongly cor- 
related systems. This is the main purpose of this pa- 
per. It has long been known that generally incommen- 
surate inhomogeneous phases arise in classical systems 
with competing short range attractive interactions and 
long range repulsive interactions. In such systems, the 
short range attractive interactions (whose physical ori- 
gin depend on the system in question) favor spatially 
inhomogeneous phases, i.e. phase separation, which is 
frustrated by long range (typically Coulomb) repulsive 
interactions. Coulomb-frustrated phase separation has 
been proposed as a mechanism for stripe phases in doped 
Mott insulators ^ 17 ' 18 and in low density electron gases ~ 
Similar ideas were also proposed to explain the struc- 
ture of the crust of neutron stars, lightly doped with 
protons ; 20 ' 21 and in soft condensed matter (e.g. block 
copolymers)^ 

In this work we will pursue a different approach and 
consider mechanisms of quantum stabilization (i.e. quan- 
tum order by disorder) of stripe-like phases in frus- 
trated quantum systems. We will specifically consider 
frustrated versions of two-dimensional quantum dimer 
models^ which provide a qualitative description of the 
physics of quantum frustrated magnets in their spin- 
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disordered phases. The phases that we will discuss here 
are essentially valence bond crystals with varying de- 
gree of commensurability and become asymptotically in- 
commensurate. Since these systems are charge-neutral, 
there are no long range interactions. As we will see 
below, quantum fluctuations resolve the high degenera- 
cies of their naive classical limit leading to a non-trivial 
phase diagram with phases with different degree of com- 
mensurability or tilt. The resulting phase diagram has 
the structure of an incomplete devil's staircase similar 
to that found in classical order-by-disorder systems such 
as the anisotropic next-nearest neighbor Ising (ANNNI) 
models 2 ^. In the regime in which quantum fluctua- 
tions are weak, which is where our calculations are sys- 
tematically controlled, only a small fraction of the phase 
diagram exhibits phases with non-trivial modulations. 
In this "classical" regime the observation of non-trivial 
phases requires fine tuning of the coupling constants. (In 
contrast, in systems with long range interactions no such 
fine tuning is needed at the classical leveljil) However, as 
the quantum fluctuations grow, the fraction of the phase 
diagram occupied by these non-trivial phases becomes 
larger. Thus, at finite values of the coupling constants, 
where our estimates are not accurate, no fine tuning is 
needed. 

Considerable progress has been made towards under- 
standing theoretically the liquid phases, including the 
proper field theoretic descriptio n 26 ' 27 which also allows 
for an analysis of the related valence bond solid phases 
with varying degrees of commensurability . 28 ' 29 Comple- 
mentary to this effort is the dynamical question of how, 
or even if, a phase with exotic non-local properties can 
arise in a system where the interactions are purely lo- 
cal. We emphasize the requirement of locality because 
many of these exotic structures have been proposed for 
experimental systems believed to be described by a local 
Hamiltonian (i. e. of the Hubbard or Heiscnbcrg type). 
An additional question is whether the exotic physics 
can be realized in an isotropic model or if it is neces- 
sary for the Hamiltonian to explicitly break symmetries. 
For some of these structures, the dynamical questions 
have been partially settled by the discovery of model 
Hamiltonian a 27 ' 30 ' 31 which stabilize the exotic phase over 
a portion of their quantum phase diagrams. While many 
of these models do not (currently) have experimental re- 
alizations, their value, in addition to providing proofs 
of principle, lies in the identification of physical mecha- 
nisms, which often have validity beyond the specific case 
considered. For example, the existence of short-range 
RVB phases was first demonstrated analytically in quan- 
tum dimer model a 23 ' 30 , the essential ingredients being ge- 
ometric frustration and ring interactions. The possibility 
of such phases being realized in spin systems was sub- 
sequently demonstrated by the construction of a number 
of spin Hamiltonians3 2 ^4' 3 £'26'37'2 8 ., including SU (2) in- 
variant model a 34 ' 35 , all of which reduce to dimcr-likc 
models at low energies. The existence of commensurate 
valence bond solid phases, with low order commensura- 



bility, in doped quantum dimer models 3 ^ has been stud- 
ied recentl y 40 ' 41 , as well as modulated phases in doped 
quantum dimer models within mean field theory- 
It is in this context that we ask the following ques- 
tion: is it possible to realize high order striped phases 
in a strongly-correlated quantum system with only local 
interactions and no explicitly broken symmetries? While 
the term stripe has been used in reference to a number 
of spatially inhomogeous states, here we use the term to 
denote a domain wall between two uniform regions. The 
presence of domain walls means that rotational symmetry 
has been broken. Fig.[T]gives examples of striped phases. 
Theories based on the formation of striped phases have 
been proposed in a number of experimental contexts, no- 
tably the high T c cuprates ^ 42 ' 43 ' 44 , where the "stripes" 
are lines of doped holes separating antiferromagnetic do- 
mains (see Fig. [IJ,). In the simplest striped states, the 
domain walls are periodically spaced (Fig.[TjD) or are part 
of a repeating unit cell (Fig. [TJ:). We use the term "high 
order striped phase" in the case where the periodicity is 
large compared to the other characteristic lengths in the 
model. 
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FIG. 1: Examples of striped phases, (a) shows stripes of 
holes separating antiferromagnetic domains. This structure 
appears in some theories of the high T c cuprates. (b) shows 
periodically spaced domain walls separating regions where the 
order parameter takes the uniform value (j>i or cf>2- (c) is 
another example where the periodicity is associated with a 
repeating unit cell. If the repeat distance becomes infinite, 
then the state is said to be incommensurate with the lattice. 
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Our central result is a positive answer to the question 
posed in the preceding paragraph. We do this by con- 
structing a two-dimensional quantum dimer model, with 
only short range interactions and no explicitly broken 
symmetries, that shows an infinite number of periodic 
striped phases in its T — phase diagram. The collec- 
tion of states forms an incomplete devil's staircase. The 
phases are separated by first order transitions and the 
spacing between stripes becomes arbitrarily large as the 
staircase is traversed. 

Before giving details of the construction, we reiter- 
ate that we are searching for (high-order) striped phases 
in a Hamiltonian with only local interactions and with- 
out explicitly breaking any symmetries. As alluded to 
earlier, a number of experimental systems where stripe- 
based theories have been proposed are widely believed 
to be described by Hamiltonians that meet these restric- 
tions. A notable example is the Emery mode l 45 ' 46 of the 
high-T c cuprates, which is a generalization of the Hub- 
bard model that includes both Cu and O sites. Since 
it is not a priori obvious that a nontrivial global or- 
dering such as a high order striped phase will/can un- 
ambiguously arise from such local, symmetric strongly 
correlated models, we include these phases in the list of 
"exotic" structures. However, in the absence of these 
restrictions, the occurrence of stripe-like phases is rela- 
tively common. Relaxing the requirement of high peri- 
odicity, we note that low order striped quantum phases 
can occur in the Bosc-Hubbard model at fractional fill- 
ings if appropriate ncxt-ncarcst-ncighbor interactions are 
added 4 ^. Relaxing the requirement of locality, we note 
that stripe phases arise naturally in systems with long 
range Coulomb interactions 1 ^ 2 -. More generally, if the 
Hamiltonian includes a term that is effectively a chemical 
potential for domain walls, and if there is a long-range 
repulsive interaction between domain walls, then we may 
generically expect striped phases where the spacing be- 
tween domain walls is large (compared with the lattice 
spacing) . 

A guiding principle in constructing models with exotic 
phases is frustration, or the inability of a system to si- 
multaneously optimize all of its local interactions. Quan- 
tum dimer models are relatively simple models that con- 
tain the basic physics of quantum frustration and have 
proven to be a useful place to look in the search for ex- 
otic phases^. This is one reason why we choose to work 
in the dimer Hilbert space. A second reason is related 
to the observation that each dimer covering may be as- 
signed a winding number and this divides the Hilbert 
space into topological sectors that are not coupled by lo- 
cal dynamics. The ground state wavefunction of a dimer 
Hamiltonian will typically live in one of these sectors (ig- 
noring the complications of degeneracy for the moment). 
As parameters in the Hamiltonian are varied, it is possi- 
ble that at some critical value, the sector containing the 
ground state will change. Such a scenario, which occurs 
even in the simplest dimer model formed by Rokhsar and 
Kivelson, is an example of a quantum tilting transition 



between a "flat" and "tilted" state. In Refs. H and H, 
it was shown, by field theoretic arguments, that when 
such a transition is perturbed, states of "intermediate 
tilt" (this will be made more precise below) may be sta- 
bilized. As we will show, these intermediate tilt states 
may be viewed as stripe-like states, of the form we are 
interested in. Taking inspiration from these ideas, our 
construction involves perturbing about a tilting transi- 
tion in a specially constructed dimer model. 

In classical Ising systems, the competition between 
nearest-neighbor and next-nearest-neighbor interactions 
is a well-known mechanism for generating incommensu- 
rate phases^ 4 -. Classical models of striped phases are 
based on two complementary principles: soliton forma- 
tion and competing interactions^ 5 -. Our quantum con- 
struction is based on analogies with two classical mod- 
els that are representatives of these two aspects: the 
Pokrovsky-Talapov model of fluctuating domain walls 4 ^ 
and the ANNNI model^l In section [HI we review the 
relevant features of these models. In section Infl we re- 
view the salient features about dimer models and tilting 
transitions. In section IIV( we give an overview of the 
construction and technical details are presented in sec- 
tion [V] and the appendices. In section I VII we discuss 
how these ideas connect to spin models, thus extending 
our "proof of principle" to systems with physical degrees 
of freedom. We discuss implications of these results in 
section fVIIl In two appendices we give technical details 
of our calculations. 



II. CLASSICAL MODELS 

Our approach builds on principles underlying stripe 
formation in classical models, where the problem is 
also referred to as a commensurate-incommensurate 
transition 2 ^. The classical models are based on two com- 
plementary principles: the insertion of domain walls and 
competing local interactions. 

A toy model relevant to the present work is the pic- 
ture of fluctuating domain walls in two dimensions, in- 
troduced by Pokrovsky and Talapov 4 ^. The walls are 
allowed to fluctuate though the ends are fixed (Fig. [2]) , 
which precludes bubbles. The free energy minimization 
is a competition between the creation energy of having 
walls, the elastic energy of deviating from the flat state, 
and the entropic benefit of allowing the walls to fluctuate. 
This theory predicts a transition from a uniform phase 
to a striped phase. The spacing between walls depends 
on the parameters of the theory (including temperature) 
and can be large compared with other length scales. 

A second model relevant to the present work is the clas- 
sical anisotropic next-nearest-neighbor Ising (ANNNI) 
model in three (and higher) dimension o 24 ' 49 . This model 
(Fig. [3J describes Ising spins on a cubic lattice with ferro- 
magnetic nearest-neighbor interactions J\ < and anti- 
ferromagnetic next-nearest-neighbor interactions J 2 > 
along one of the lattice directions. A key feature of the 
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FIG. 2: (a) ground and (b) excited states of the Pokrovsky- 
Talapov model of fluctuating classical domain walls. This is a 
two-dimensional anisotropic model where domain walls form 
along the y direction and separate regions where the order 
parameter is uniform. While the domain walls cost energy, 
they are allowed to fluctuate, which carries entropy. For a 
range of parameters, the domain walls are actually favored 
by the free energy minimization. 




J1 

FIG. 3: The anisotropic next nearest neighbor Ising (ANNNI) 
model. Ising spins lie on the points of a d-dimensional cube. 
Nearest neighbor interactions are ferromagnetic (Ji < 0. 
Along one of the directions, we also have antiferromagnetic 
next nearest neighbor interactions J2 > 0. 



ANNNI Hamiltonian is a special point J1/J2 = 2 where 
a large number of stripe-like states are degenerate at zero 
temperature. As the temperature is raised, the compe- 
tition between J\ and J2 causes an infinite number of 
modulated phases to emerge from this degenerate point. 
The phase diagram in the low T limit was studied analyt- 
ically in Ref. [2J using a novel perturbative scheme where 
the existence of higher order phases was established at 



successively high orders in the perturbation theory. Nu- 
merical studies at higher temperatures^ indicated that 
incommensurate phases occur near the phase boundaries. 
Therefore, the collection of phases form an incomplete 
devil's staircas o 25 i 50 . A quantum version of the ANNNI 
model was studied in Rcfs. [5l| andfU 

The phase diagram of our model is similar to that of 
the ANNNI model and our analytical methods are similar 
in spirit to Ref. [24|. However, the basic physics of our 
model corresponds more clearly to a quantum version 
of the energy-entropy balance occuring in the fluctuating 
domain wall picture. We now discuss one more ingredient 
of the construction before putting the pieces together in 
section [IV] 



III. QUANTUM DIMER MODELS AND 
TILTING TRANSITIONS 

A hard-core dimer covering of a lattice is a mapping 
where each site of the lattice forms a bond with exactly 
one of its nearest neighbors. Each dimer covering is a 
basis vector in a dimer Hilbert space and the inner prod- 
uct is such that different dimer coverings are orthogonal. 
Quantum dimer models are defined on this dimer Hilbert 
space through operators that manipulate these dimers in 
ways that preserve the hard-core condition. These mod- 
els were first proposed as effective descriptions of the 
stron g couplin g regime of quantum spin systems^ and 
Refs.iHiii discuss ways in which this correspondence 
can be made precise. 

The space of dimer coverings can be subdivided into 
different topological sectors labelled by the pair of wind- 
ing numbers (W x ,W y ) defined in Fig. [4] The winding- 
number is a global property in that it is not affected 
by any local rearrangement of dimers. In particular, for 
any local Hamiltonian, the matrix element between dimer 
coverings in different sectors will be zero. 

For a given local Hamiltonian, the ground state of the 
system will typically lie in one of the topological sec- 
tors. A common occurrence in dimer models is a quan- 
tum phase transition in which the topological sector con- 
taining the ground state changes. Such an occurrence is 
called a tilting transition because a dimer covering of a 2d 
bipartite lattice may be viewed as the surface of a three- 
dimensional crystal through the height representation^. 
In this language, the different topological sectors cor- 
respond to different values for the (global) average tilt 
of the surface. The correspondence between dimers and 
heights is reviewed in Fig. [5] but for the present purpose, 
it is sufficient to define the "tilt" of a dimer covering 
as its "winding number per unit length" . The simplest 
dimer model introduced in Ref. [23l . has a tilting tran- 
sition between a flat state (zero tilt) and the staggered 
state, which is maximally tilted (Fig. [SJd). At the criti- 
cal point, called the Rokhsar-Kivelson or RK point, the 
Hamiltonian has a ground state degeneracy where all tilts 
are equally favored. 
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(a) 



FIG. 4: Winding numbers: draw a reference line that extends 
through and around (due to the periodic boundary condition) 
the system and label the vertical lines of the lattice alternately 
as A and B lines. For any dimer configuration, we may de- 
fine, with regard to this reference line, the winding number 
W x = Na — Nb, where Na is the number of A dimers inter- 
secting the line and similarly for Nb ■ We can similarly draw a 
vertical line and define a similar quantity W y . Note that this 
particular construction works for a 2d bipartite lattice. For 2d 
non-bipartite lattices, the construction is simpler: count the 
total number of dimers intersecting the horizontal and verti- 
cal reference lines and there are four sectors corresponding to 
whether W x ,y is even or odd. 



The recognition of the Rokhsar-Kivelson dimer 
model as a tilting transition led to a field theoretic 
descriptio n 26 ! 54 of the RK point based on a coarse- 
grained^ 5 - version of the height field (Fig. [5]). The sta- 
bility of this field theory was studied in Refs. andl29l. 
These studies showed that by tuning a small perturba- 
tion and non-perturbatively adding irrelevant operators, 
it is possible to make the tilt vary continuously from a 
flat state to the maximally tilted state. In addition, it 
was observed that the system has a tendency to "lock- 
in" at values of the tilt commensurate with the under- 
lying microscopic lattice, the specific values depending 
on details of the perturbation. It was also noted that 
while a generic perturbation would make the transition 
first order—, for a sufficiently small perturbation, the 
correlation length was extremely large 2 ^ which, it was 
argued, justified the field theory approach nonetheless. 
Therefore, the generic effect of perturbations would be 
to smoothen the Rokhsar-Kivelson tilting transition by 
making the system pass through an incomplete devil's 
staircase of intermediately tilted states. One may sus- 
pect that the structure of the field theory, including the 
predictions of Refs. [28| andU^, would hold for a broader 
class of tilting transitions. In particular, one may con- 
sider the case where the critical point is merely a point 
of large degeneracy where all tilts are favore d 56 ! 57 , which 
is analogous to the classical ANNNI model. 

The relevance of all of this to the present work is most 
easily seen by Fig. [9l which shows the simplest examples 
of states that have intermediate tilt (i. e. winding num- 
ber). These are stripe-like states having a finite density 
of staggered domain walls and more general tilted states 
may be obtained by locally rearranging the dimers. The 
preceding discussion suggests that these kinds of struc- 
tures arise naturally when quantum dimer models are 
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FIG. 5: Sample dimer configurations with corresponding 
height mappings. The height mapping involves assigning in- 
tegers to the squares of the lattice in the following manner. 
Divide the bipartite lattice into A and B sublattices. Assign 
zero height to a reference square and then moving clockwise 
around an A site, the height increases by one if a dimer is 
not crossed and decreases by three if a dimer is crossed. The 
same rule applies moving counterclockwise about a B site. 
The integers correspond to local heights of a crystal whose 
base lies on the page. In these examples, the lower square 
in the second column is taken as the reference square, (a) 
Dimers are arranged in columns corresponding to a surface 
that is flat on average (though there are fluctuations at the 
lattice scale) . (b) Dimers are staggered and the corresponding 
surface is maximally tilted, (c) There average tilt is nonzero 
due to the staggered strip in between the flat columnar re- 
gions, (d) Going from left to right, the surface initially falls 
and then rises giving an average tilt of zero. This is because 
the two staggered regions have opposite orientation. 



perturbed around a tilting transition. This observation 
will guide the construction outlined in the next section. 



IV. OVERVIEW OF STRATEGY 

We now combine various ideas presented in sections |TT] 
and Mil to construct the promised quantum model. In 
the present section, we present an overview of the con- 
struction with details and subtleties relegated to section 
fVl The basic idea is to construct a quantum dimer model 
with a tilting transition and then to appropriately per- 
turb this model to realize a staircase of striped states. 
We design the unperturbed system to have a large degen- 
eracy at the critical point, with each of the degenerate 
states having a domain-wall structure. The perturba- 
tion will effectively make these domain walls fluctuate 
and the degeneracy will be lifted in a quantum analog of 
the energy-entropy competition that drives the classical 
Pokrovsky-Talapov transition. Using standard quantum 
mechanical perturbation theory, we will obtain a phase 
diagram similar to the classical 3D ANNNI model and 
will find that phases with increasingly long periodicities 
will be stabilized at higher orders in the perturbation 
theory. This mathematical approach is similar in spirit 
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to the analysis of the classical ANNNI model in Ref. l24l . 

A simple, rotationally invariant Hamiltonian with a 
tilting transition is given by: 

H o = 2-a(lB)(Bl + l0XBI + lm)(m| + im)(ml) + 
2-b(|nXn|+|nXn|) 

(4.1) 

This model displays a first-order transition between a 
columnar and fully staggered state at a very degenerate 
point, where 2a = b. In principle, we may perturb this 
model with an off-diagonal resonance term and expect 
phases with intermediate tilt (and possibly other exotic 
phases) on the general grounds discussed in the previous 
section. However, it is difficult to make precise state- 
ments about the phase diagram even for such fairly sim- 
ple models. We will study a slightly constrained version 
of this model that is convenient for making analytical 
progress. 

We construct the quantum dimer model in two steps. 
First, we construct a diagonal parent Hamiltonian Hq 
(Eq. 15. ip where the ground states are separated from ex- 
cited states by a tunably large gap. Hq will not break any 
lattice symmetries, but the preferred ground states will 
spontaneously break translational and rotational symme- 
try. In particular, we design Hq to select ground states 
having the domain wall structure shown in Fig. [BJ In 
these states, the dimers arrange themselves into stag- 
gered domains of unit width separating columnar regions 
of arbitrary width. The columnar dimers are horizontal if 
the staggered dimers are vertical (and vice versa) and the 
staggered strips come in two orientations. Notice that 
the fully columnar state (a columnar region of infinite 
width) is included this collection but the fully staggered 
state (Fig. (Sb), which appears in the Rokhsar-Kivelson 
phase diagram, is not. In the following, we will typically 
draw the staggered strips as vertical domain walls but 
the horizontal configurations are equally possible. 




FIG. 6: A typical domain wall state selected by Ho- These 
states break translational and rotational symmetry. The stag- 
gered strips, which are one column wide and may have one 
of two orientations, are like domain walls separating colum- 
nar regions, which may have arbitrary width. When a — b in 
Eq. 15.21 the set of these states spans the degenerate ground 
state manifold of Hq. 



In analogy with the ANNNI model, Hq is designed so 
that all of these domain wall states are degenerate when 
the couplings arc tuned to a special point. Away from 
this point, the system will enter either a flat or a tilted 
phase. The unperturbed phase diagram is sketched in 
Fig.m 




columnar ; herringbone 



(highly degenerate) 




FIG. 7: Ground state phase diagram of the parent Hamilto- 
nian Ho as a function of the parameter a — b (see Eq. I5.1[l . In 
these states, the dimers may only have two attractive bonds. 
When a — b — 0, the states of Fig. [5] are degenerate ground 
states. Away from this point, the system enters a state where 
dimers either maximize or minimize the number of staggered 
interactions. The maximally staggered configuration is com- 
monly called the "herringbone" state. 



The second step of the construction is to perturbativcly 
add a small, non-diagonal, resonance term tV: 

-tv = 2-t(|0 xgi+|ra>(m |+h.c.) 

(4.2) 

The sum is over all plaquettes in the lattice. Depend- 
ing on the local dimer configuration of the wavefunction, 
the individual terms in this sum will cither annihilate the 
state or flip the local cluster of dimers as shown in Fig. [8] 
The action of this operator on the domain wall states 
(Fig. is confined to the boundaries between staggered 
and columnar regions and effectively makes the domain 
walls fluctuate. Notice that Eq. 14.21 is equivalent to two 
actions of the familiar Rokhsar-Kivelson two-dimer reso- 
nance term. We expect the basic conceptual argument to 
apply for a more general class of perturbations, includ- 
ing the two-dimer resonance, but we consider the specific 
form of Eq. 14.21 to simplify certain technical aspects of 
the calculation. We will elaborate on this more in the 
next section. 

The degenerate point of Fig.[7]may be viewed as the de- 
generacy of an individual vertical column having a stag- 
gered or columnar dimer arrangement. The perturba- 
tion 14.21 lifts this degeneracy by lowering the energy of 
configurations with domain walls by an amount of order 
~ Lt 2 per domain wall, where L is the linear size of the 
system. Therefore, the system favors one of the states 
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FIG. 8: One of the terms in the operator of Eg. 14.21 will flip 
the circled cluster as shown. 




(e) (f) 

FIG. 9: Examples of ideal tilted states. In these states, the 
domain walls have the same orientations and are uniformly 
spaced. The notation [In] denotes the state where one stag- 
gered strip is followed by n columnar strips and so on. It 
is understood that [In] collectively refers to the above states 
and those related by translational, rotational, and reflection 
(i.e. switching the orientation of the staggering) symmetries. 
The examples drawn here, where it is understood that what 
we are seeing is part of a larger lattice, are (a) [11] (b) [12] 
(c) [13] (d) [14] (e) [15] (f) [16]. 



with a maximal number of domain walls and for techni- 
cal reasons discussed below, will choose the one having 
maximal tilt: the [11] state in Fig. [9^,. However, the de- 
generacy between columnar and staggered strips will be 
restored by detuning H from the t = degenerate point. 
This implies the second order phase diagram sketched in 

Fig. ma 

Fig. [10] is correct up to error terms of order t 4 . To 
this approximation, states having tilt in between the [11] 
state and the (flat) columnar state are degenerate on the 
phase boundary. Physically, this boundary occurs when 
the energy from second order processes which stabilize 
the staggered domains is precisely balanced by the energy 
of a columnar strip. This degeneracy will be partially 
lifted by going to higher orders in perturbation theory. 
We find that at fourth order, a new phase is stabilized 
in a region of width ~ t 4 between the columnar and [11] 
phases. This new phase is the one which maximizes the 




FIG. 10: (Color online) Ground state phase diagram of H = 
Ho + tV from second order perturbation theory. 



number of fourth order resonances shown in Fig. [TTJ and 
is the [12] state (Fig. (9Jd). The corrected phase diagram 
is given in Fig. 1121 



initial state 




FIG. 11: (Color online) The excited state on the right is ob- 
tained from the initial state by acting twice with the pertur- 
bation in Eq. 14.21 




FIG. 12: (Color online) Ground state phase diagram of H = 
Ho + tV from fourth order perturbation theory. The width of 
the [12] phase is order t 4 . 

Fig. H2l is accurate up to corrections of order t 6 . To 
this approximation, on the boundaries, states with tilts 
in between the bordering phases are degenerate. These 
degeneracies will be partially lifted at higher orders in 
the perturbation. At higher orders, there will be reso- 
nances corresponding to increasingly complicated fluctu- 
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FIG. 13: (Color online) Ground state phase diagram of H = 
Ho + tV from nth order perturbation theory. The width of 
the [In] phase is of order t 2n . 




FIG. 14: (Color online) The boundaries of the [In] sequence 
will typically open into finer phases and subsequently the fine 
boundaries can themselves open. While the detailed structure 
depends on parameters in the model, generically we expect an 
incomplete devil's staircase to be realized. In the figure,we 
have explicitly drawn the opening of the[ll]-[12] boundary 
but the other boundaries will behave similarly. 



ations of the staggered lines but at nth. order, the com- 
petition between the [l,n-l] and columnar phases will 
always be decided by the nth order generalization of the 
long resonance in Fig. 111! The competition will stabilize 
a new [In] phase in a tiny region of width ~ t 2n between 
the [l,n-l] and columnar phases resulting in the phase 
diagram of Fig. [TH 

We will also find that at higher orders, the individ- 
ual boundaries of the [In] sequence will themselves open 
revealing finer phase boundaries, which themselves can 
open. This leads to the generic phase diagram sketched in 
Fig- [HI The steps in the [In] sequence that are stabilized 
depend on the values of parameters in the Hamiltonian. 
However, the conclusion of arbitrarily long periods being 
realized is robust. The fine structure of how the [l,n-l]- 
[ln] boundaries open is less certain because the depen- 
dence on parameters is more intricate and increasingly 
complicated resonances need to be accounted for. How- 
ever, general arguments indicate that the boundaries will 



open and even periods incommensurate with the lattice 5 ^ 
are likely to occur in the model, though such states will 
not be seen at any finite order of our perturbation theory. 

In the terminology of commensurate-incommensurate 
phase transitions, the [In] sequence forms a (harmless) 
staircase with a "devil's top-step" 24 ' 25 . With the open- 
ings of these boundaries, beginning in the [11] state and 
moving left in Fig. [TJ] for t ^ 0, the system traverses an 
incomplete devil's staircase of periodic states. The sub- 
sequent steps in the staircase have progressively smaller 
tilts culminating in the flat columnar state. The phase 
boundaries are first order transitions. This phase dia- 
gram is similar to what is seen in the classical ANNNI 
model, where the transitions are driven by thermal fluc- 
tuations. 

We make two remarks before launching into the cal- 
culation. First, when we refer to the "[11] phase" (for 
example) what we precisely mean is that in this region, 
the ground state wavefunction is a superposition of dimcr 
coverings that has relatively large overlap with the state 
in Fig. [5^ and much smaller overlaps (of order t 2 , t 4 etc. 
) with excited states obtained by acting on Fig. [9^ with 
Eq. 14.21 The coefficients follow from perturbation the- 
ory. Second, since Fig. [T3] is obtained using perturba- 
tion theory, we can be confident that this describes our 
system only in the limit where t is small. In the classi- 
cal ANNNI model, numerical evidence indicates that as 
the small parameter (the temperature in that case) is in- 
creased, the phase boundaries close into Arnold tongue 
structures. We do not currently know if this will occur 
in our model as t increases. 



V. DETAILS 

In this section, we construct a Hamiltonian using the 
strategy outlined in the previous sections. The Hamilto- 
nian H = Hq + tV consists of a diagonal term Ho and 
an off-diagonal term tV, which we treat perturbatively 
in the small parameter t. 



A. Parent Hamiltonian 

Our parent Hamiltonian Hq is the following operator: 

q, = 2-a(6XBMBXSMm><mMmXml) + 
2-b(|P)(P | +!□><□ |) + 

2c(|EBXE5l +3more) + 

2d(iEBXEEl +7m ° re ) + 
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(a) 






(d) 



FIG. 15: (Color online) (a) and (b) are the two ways in which 
a dimer can participate in three attractive bonds, (c) is the 
one way in which a dimer can participate in four attractive 
bonds. The attractive bonds are shown by the blue (dark 
gray) arrows. However, these configurations also involve re- 
pulsive interactions, which are shown in red (light gray) , from 
terms c and d in the Hamiltonian. (d) is an example of a state 
where every where every dimer has only two attractive bonds 
but with some dimers the two bonds are of different types. 
These "kinks" in the staggered domain walls involve an en- 
ergy cost from term d as indicated by the red (light gray) 
arrows. 



2p(| mi )(□=□ | + 3 more) 

2q(lFEBxFEBl +7more ) 



(5.1) 

The coefficients a, b, c, d, p, and q are positive numbers. 
The symbols used in Eq. 15. H are projection operators re- 
ferring to configurations of dimers on clusters of plaque- 
ttes and the sums are over all such clusters. The notation 
"3 more" , etc. refers to the given term as well terms re- 
lated to it by rotational and/or reflection symmetry; in 
terms a and b these other terms are explicitly written. 
Notice that this Hamiltonian is a sum of local operators 
and does not break any symmetries of the underlying 
square lattice. 

Terms a and b are attractive interactions favoring stag- 
gered and columnar dimer arrangements respectively and 
we study Eq. 15. II near a = b. Terms c and d are repulsive 
interactions and if c, d > a, b, the dimers prefer domain 
wall patterns (Fig. [6]). Terms p and q are repulsive in- 
teractions which determine the phases on the staircase. 
If these terms are sufficiently large^ compared to a, 6, c 
and d, the staircase will include phases with arbitrarily 



long periods. 

We begin by showing that when a = b, the ground 
states of Hq are the domain wall states of Fig. [H We do 
this by showing that competitive states must have higher 
energy. In the domain wall states, every dimer partici- 
pates in exactly two attractive interactions and no repul- 
sive interactions. The only way to achieve a lower energy 
is for some dimers to participate in three or four attrac- 
tive interactions. This involves local dimer patterns of 
the form shown in Fig. [T5bbc. In Fig. H5b . the central 
dimer participates in one columnar and two staggered in- 
teractions but also two repulsive interactions from term 
d in Eq. 15.11 Similarly, Figs. [T5b and [T5b show that if 
a dimer participates in more than two staggered inter- 
actions, the extra bonds are penalized by term c. If we 
require c,d> a, 6, these patterns will result in higher en- 
ergy states as the repulsive terms nullify the advantage of 
having extra attractive bonds. This also explains why in 
Fig. O the staggered strips have unit width and the stag- 
gered and columnar dimers have opposite orientation. 

Of the states where every dimer has two attractive 
bonds, the states where some dimers have two bonds 
of different type will also have higher energy as shown 
in Fig. [T5H . Of the remaining states, it is readily seen 
that states where every dimer participates in either two 
a bonds or two b bonds, and where there are some b 
bonds, must be of the domain wall form. The only other 
possibility is the "herringbone state" where every dimer 
has two a bonds (see Fig. [7|). The latter states are part 
of the degenerate manifold at a — b but are dynamically 
inert because in this state it is not possible to locally ma- 
nipulate the dimers (without violating the hard-core con- 
straint). This establishes that when a = b, the ground 
states have the domain wall form. It is also clear that 
when a < b, the system will maximize the number of b 
bonds and when a > b, the number of a bonds. There- 
fore, we obtain the zero temperature phase diagram in 
Fig. [7] 

It is useful to see this formally by calculating the en- 
ergy of each domain wall state. For concreteness, we 
assume the translational symmetry is broken in the x di- 
rection. A configuration with N s staggered strips and N c 
columnar strips has energy: 



E(N S ,N C ) 



-aN y N s - bN y N c 



+ (b- a)N y N s (5.2) 



where N x and N y are the dimensions of the lattice (the 
lattice spacing is set to unity) . We have used the relation 
N s + N c = ^ . When a = b, the domain wall states are 
degenerate and the energy scales with the total number 
of plaquettes N x N y . If a < 6, the system prefers the min- 
imal number of staggered strips, which is the columnar 
state. If a > 6, the herringbone configuration has lower 
energy than any domain wall state. Notice that all of 
these states are separated by a gap of order a or & from 
the nearest excited states obtainable by local manipula- 
tions of dimers. Since we will be interested mainly in the 



difference a — b, the individual size of a (or b) , which sets 
the scale of this gap, can be made arbitrarily large. 

B. Perturbation 

We now consider the effect of perturbing the parent 
Hamiltonian (|5.ip with the non-diagonal resonance term 
given in Eq. 14.21 

-tv = s-t(|0 >(Bl+|ra><m |+h.c.) 

(5.3) 

We assume t << l,a,b and consider f as a small pa- 
rameter in perturbation theory. We examine how the 
degeneracies in the t = phase diagram (Fig. [7]) get 
lifted when t ^ 0. The technical complications of degen- 
erate perturbation theory do not arise because different 
domain wall states are not connected by a finite number 
of applications of this operator. Eq. 15.31 is equivalent to 
two applications of the familiar two-dimer resonance of 
Rokhsar-Kivclson. The mechanism we now discuss can, 
in principle, be made to work for even this two-dimer 
term, but there arc additional subtleties which will be 
mentioned. 



1. Second order 

An even number of applications of operator 15.31 arc 
required to connect a domain wall state back to itself. 
Therefore, to linear order in t, the energies of these states 
are unchanged. To second order in t, the energy shift of 
state \n) is given by: 

E n = €n _ t 2J2 V ™V mn +Q(t4) (54) 

e n is the unperturbed energy of state \n) as given by 
Eq. 15.21 The primed summation is over all dimer cov- 
erings except the original state \n). The terms in the 
sum which give nonzero contribution correspond to states 
connected to the initial state by a single flipped cluster. 
These terms may be interpreted as virtual processes tak- 
ing the initial state to and from higher energy intermedi- 
ate states, which may be viewed as quantum fluctuations 
of the staggered lines. 

The resonance energy of a particular staggered line de- 
pends on how the dimers next to the line are arranged. 
Fig. 1161 shows the three possibilities that give denomina- 
tors that may enter Eq. 15.41 The staggered line may be 
of type (1) next to a columnar region of greater than 
unit width (Fig. H6ll). or separated by one columnar 
strip from another staggered line of type (2), having the 
same orientation (Fig. [T61-2) or type (3), having opposite 
(Fig.[TU-3) orientation. From the figure, we sec that cases 
(2) and (3) involve higher energy intermediate states than 
case (1) but allow for a denser packing of lines. 
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FIG. 16: (Color online) The three types of intermediate states 
obtained by acting once with the perturbation 15.31 on the do- 
main wall states. The blue/dark gray (red/light gray) arrows 
denote attractive (repulsive) interactions that are present in 
the initial (left) or excited (right) state but not both. The 
circled cluster of the excited state of (3) is another such in- 
teraction. (1) is a staggered line next to a columnar region 
of greater than unit width. (2) and (3) are staggered lines 
separated by one columnar strip from another staggered line 
of the same and opposite orientation respectively. Notice that 
relative to the excited state of (1), the excited state of (2) has 
2 additional attractive a bonds and 2 additional repulsive c 
bonds. Similarly, excited state (3) has 1 additional attractive 
b bond, 2 additional repulsive d bonds, and a repulsive p in- 
teraction. Since c,d > a, b, cases (2) and (3) involve higher 
energy intermediate states than (1). 



If p is sufficiently large, states with type (3) lines will 
be disfavored as ground states. Ignoring such states, we 
update Eq. 15.21 to include second order corrections. For 
convenience, we set c = a which removes the distinction 
between cases (1) and (2), 

N N t 2 
E(N S ) = -l^' +{b -a- w ^)N v N. 

(5.5) 

We may use this to update the ground state phase dia- 

.2 

gram. If b > a + 4d+2b , the system is optimized when 

N s = 0, which is the columnar phase. If b < a + 4( / +2b , 
the best domain wall state is the one with maximal stag- 
gering but without case (3) lines. This corresponds to 
the [11] state (Fig. [5]a,) in which every staggered strip 
has the same orientation. As a — b is further increased, 
the herringbone state will eventually be favored. The 
boundary between the [11] and herringbone state may 
be determined by comparing energies. The energy of the 
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A2 N 



[11] 



columnar 



herringbone 



b-a 



/ 

highly degenerate 



FIG. 17: (Color online) Second order phase diagram when 
c > a. The [11] phase has width ~ t 2 and the A2 phase 
has width ~ (c — a)t 2 . On the [llJ-^2 and the ^-columnar 
boundaries, intermediate domain wall states are degenerate 
as described in the text. 



[11] state is: 

E {11] = (-b 



t- 



Ad+2b 



N X Ny 



+ 0(t 4 ) (5.6) 



while the herringbone state has energy Eh = 
From this, it follows that the [11] state will be favored 
when: 



4d+2b 



< b < a 



t 2 



Ad+2b 



(5.7) 



while the herringbone state will occur when b < a— 4d * +2b . 

Therefore, up to corrections of order t 4 , the system 
has the phase diagram shown in Fig. [101 Because the 
coefficient of N s in Eg. 15.51 is zero on the [ll]-columnar 
boundary, we have that the [11] state, columnar state, 
and any domain wall state with intermediate tilt (that 
contains only lines of type (1) and (2)) are degenerate 
on the boundary. In contrast, on the [ll]-herringbonc 
boundary, only the two states are degenerate. 

In Appendix E] the more general case, where c > a, 
is discussed and the resulting phase diagram is shown in 
Fig. [T7] An additional phase is stabilized in a region of 
width ~ (c — a)t 2 (assuming |c — a\ is finite) between 
the [11] and columnar states. In this new phase, labelled 
A 2 , any state where adjacent staggered lines are sepa- 
rated by two columns, including the [12] state (Fig. [3b), 
is a ground state. These are the states which maximize 
the number of type (1) staggered lines (Fig. [TBT-l). On 
the [11]-j42 boundary, intermediate states where adjacent 
staggered lines are separated by one or two columns (and 
where there are no type (3) lines) are degenerate. On the 
^-columnar boundary, states where adjacent staggered 
lines are separated by at least two columns are degener- 
ate. 

Fig. [17] may be understood intuitively by noting that 
at a = b and t = 0, a staggered strip has the same energy 
as a columnar strip. The resonance terms lower the ef- 
fective energy of a staggered strip and since the [11] state 



involves the most staggered strips, its energy will be low- 
ered the most. As b increases to the point where the 
degeneracy between columnar strips and type (2) lines is 
restored, the system will prefer to maximize the number 
of type (1) lines, which are individually more stable but 
loosely packed. This is the transition to the A2 phase. As 
b is increased further, the degeneracy between columnar 
strips and type (1) lines is restored and there is a tran- 
sition to the columnar state. Note that if \c — a\ is very 
large, [11]-A 2 boundary can occur in the a > b region. 

Before proceeding with the calculation, we clarify two 
points of potential confusion. First, when we say the 
"[11] state is stabilized" over part of the phase diagram, 
what we precisely mean is that the ground state wave- 
function is a superposition of dimer coverings that has 
relatively large overlap with the literal [11] state of Fig.[9K 
and much smaller overlaps (of order t ) with the excited 
states obtained by acting on the [11] state once with the 
perturbation l5.3l fFig.l8l). We use the notation [11] to de- 
note both the perturbed wavehmction, which is an eigen- 
state of the perturbed Hamiltonian and the literal [11] 
state, which is an eigenstate of the unperturbed Hamil- 
tonian. Second, the phase boundaries are based on a 
competition between a resonance term, which is a quan- 
tum version of "entropy" , and part of the zeroeth order 
piece which, continuing the classical analogy, is like an 
internal energy. This does not contradict the spirit of 
perturbation theory because the full zeroeth order term, 
_j J NyN SL _|_ ^ _ a )7\r y A s , is always larger than the second 
order correction. 

A similar phase diagram will occur (though at fourth 
order in perturbation theory) if we use the two-dimcr 
move of Rokhsar-Kivclson. However, the bookkeeping 
will be more complicated because resonances will be able 
to originate in the interior of the columnar regions in- 
stead of just at that the columnar-staggered boundaries. 
This may be compensated by tuning 6, which will merely 
move the boundaries, or by adding appropriate (local) 
repulsive terms to the parent Hamiltonian. 



2. Fourth order 



We concentrate on Fig. [17] as it is more generic than 
the fine-tuned c = a case of Fig. [TO] In cither case, we 
expect the degeneracies on the phase boundaries to be 
partially lifted by considering higher orders in perturba- 
tion theory. In this section, we focus on the ^-columnar 
boundary, where adjacent lines are separated by at least 
two columnar strips. 

We return to the perturbation series for the energy 
(Eq. I5.4[) , this time keeping terms up to fourth order in 
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the small parameter. 



m 

E 



'7X717. 



nm » 77X71 



(e m - e„)(e; - e n )(e fc - e n ) 

'Tint* run *nl *ln 



inl 



(e m - en) 2 (ei - e«) 



+ o(r) (5. 



As usual, the primes denote that the sum is over all 
states except \n) . The two fourth order terms correspond, 
in conventional Raylcigh-Schrodingcr perturbation the- 
ory, to the corrections to the energy expectation value 
(x/)\H\'ip) and wavefunction normalization We will 

refer to the corresponding terms as "energy" and "wave- 
function" terms. As in the second order case, we may 
view the terms in the fourth order sums as virtual reso- 
nances connecting the initial state to itself via a series of 
high energy intermediate states. For this reason, we will 
refer to the summands as (fourth-order) "resonances" . 

Most terms in the double sum in Eq. 15.81 correspond 
to resonances between disconnected clusters (Fig. [PSf . 
Referring to the figure, we use the term "disconnected" 
to indicate that there are no interaction terms connecting 
the dimers of clusters 1 and 2. While the number of such 
resonances scales as the square of the system size, the 
contributions from the energy and wavefunction terms in 
Eq. 15.81 precisely cancel for these disconnected clusters. 
The details of this cancellation are discussed in Appendix 

m 

The remaining fourth order resonances are extensive 
in number and may be grouped into three categories. 
In the first category, shown in Fig. 1191 are resonances 
associated with a single staggered line and the number 
of such resonances in the system is proportional to the 
number of lines. We refer to these resonances as "self- 
energy corrections" . 

In the second category, shown in Fig. [201 there are 
resonances that contribute to the effective interactions 
between adjacent lines. These resonances occur only in 
states where lines are separated by two or fewer colum- 
nar strips. Because we are interested in the ^-columnar 
boundary, the only processes to consider are the ones 
shown in the figure. The purpose of terms p and q in 
Eq. 15. II is to control the processes in Figs. [20k and b re- 
spectively. The net contribution of these resonances in- 
volve both energy and wavefunction terms. For example, 
the contribution of Fig. [20b to the energy is: 



2t l 



e = 



2(4d + 2o) ; 



1 



1 



2q— a 
2(4d+2b) 



(5.9) 



and likewise for Fig. [20k (replace 2q— a with 2p—b). If q > 
a/2, the net contribution of this resonance is repulsive. 

The third type of resonance is the long resonance, 
shown in Fig. 1211 These resonances occur in the energy 



term of Eq. 15.81 but do not have corresponding pieces 
in the wavefunction term. Therefore, these resonances 
will always lower the energy though, as the figure indi- 
cates, the precise amount depends on the way the lines 
are spaced. 

initial state (0) 
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excited state (2) 



excited state (12) 

FIG. 18: (Color online) Most of the fourth order terms in 
Eq. 1 5. 8 1 involve "disconnected" clusters of dimers. In this fig- 
ure, the perturbation connects the initial state (0) to excited 
states, labelled (1) and (2), depending on whether cluster 1 or 
2 has been flipped. Acting again with the perturbation con- 
nects to an excited state, labeled (12), where both of these 
clusters are flipped. Acting two more times with the perturba- 
tion brings us back to the initial state (0) via either of excited 
states (1) or (2). Such terms are called disconnected because 
there are no interactions (in Eq. I5,l[l between the dimers of 
clusters 1 and 2. The figure depicts a particular resonance 
from the energy term in Eq. 15.81 There is an analogous contri- 
bution from the wavefunction term which is a product of the 
second order processes (0) — » (1) — > (0) and (0) — » (2) — > (0). 
While the number of such disconnected terms scales as iV 2 , 
where N = L x L y is the system size, these resonances do not 
contribute to the energy because the contributions from the 
energy and wavefunction terms exactly cancel. 



An immediate implication of these resonances is the 
lifting of the degeneracy of the phase. All of these 
states have the same number of lines so will receive the 
same self-energy contribution (Fig. H"9j) . If we choose p, q 
large compared to a, 6, c and d, then the repulsive con- 
tribution from the interaction resonances (Fig. [20[) is es- 
sentially determined by the wavefunction term, which is 
the same for all of the A 2 states. The degeneracy is bro- 
ken by the long resonances because in the [12] state, only 
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FIG. 20: (Color online) These are examples of fourth order 
resonances which are effective interactions between lines. At 
fourth order, such interactions are possible only when lines 
are separated by two or fewer columns. On the ^-columnar 
boundary, resonances (a) and (b) are the only processes to 
consider. These involve terms p and q in the Hamiltonian, as 
indicated by the circles. 
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FIG. 19: (Color online) These are examples of fourth order 
self-energy resonances. Each resonance is confined to a single 
line and the number of resonances in the system is propor- 
tional to the number of lines. In resonances such as (a), there 
are interactions connecting dimers of two flipped clusters on 
the same line. Terms such as (b) arise from the wavefunction 
term in Eq. l5.8l but have no analogous processes in the energy 
term to cancel against. 
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FIG. 21: (Color online) These are the long fourth order res- 
onances which occur in states where lines are separated by 
two or more columnar strips. These processes are always sta- 
bilizing though the amount depends on the environment of 
the line as in the second order (see Fig. I16[) . The resonance 
in (a) is strongest because it involves the lowest energy inter- 
mediate state but (b) and (c) allow for a denser packing of 
lines. Resonance (c) is especially suppressed due to term p in 
Eq.EH 



Fig. [2Tb processes occur while in the other states, 
some of the resonances are the suppressed Fig. [21b vari- 
ety. Therefore, what was seen as merely an Ai phase at 
second order is revealed, on closer inspection, as a [12] 
phase. 

To investigate the degeneracy of what we now recog- 
nize as the [12] -columnar boundary, it is useful to update 



Eq. 15.51 to include fourth order corrections 

NyN x 



E 



+ (*-«- ^— + at*)N y N a 



- (3t*N y N sa -7t 4 N sb + 0(t 6 ) (5.10) 
where N s is the total number of staggered lines and N sa ^ 
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is the number of staggered lines having the environment 
of Fig. l2"Th(b). We ignore states with arrangements like 
Fig. l2"Tb since they are disfavored as ground states, a is 
a constant, which may be calculated but whose value is 
unimportant, containing the contribution of fourth order 
self energy terms. (3 = (4d+2b) . (2(4d 1 +2b)+4( ^ Q)) > is 
the contribution of the most favorable long fourth order 
resonances (Fig. 121b ) whose number is proportional to 

N sa ■ 7 = (4d+2by^(2(4d+2b)+4{d-a)+2(c-a)) ~ ^ e ( e given 

by Eq. 15. 9jl is proportional to N s b and includes the con- 
tributions of Figs. [20b andl2"Tb. Note that 7 < (3 and the 
sign of 7 is determined by the size of q. For convenience, 
we assume q large enough that 7 < 0. 

We may use Eq. 15.101 to correct the phase diagram. 
Similar to the second order case, as b is increased, the 
extra stability of the staggered strips in the [12] state be- 
comes eventually balanced by the zeroeth order energy of 
the columnar strips. When this occurs, the system will 
prefer a state with fewer lines that are individually more 
stable. In particular, the states we may label A3, where 
adjacent lines are separated by three columns, maximize 
the number of favorable long resonances (Fig. l2Th ) with- 
out incurring any of the repulsive fourth order penalties 
(i. e. the analog of Fig. [20] would be a disconnected res- 
onance). In the next section, higher order perturbation 
theory will show that this A 3 phase is actually a [13] 
phase (Fig. [3b) so we begin using the [13] label immedi- 
ately. As b is increased further, the system will enter the 
columnar state. 

The result is the phase diagram in Fig. [53] The phase 
boundaries are determined by comparing energies. Ig- 
noring corrections of order i 6 , we have the following: 
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Comparing these expressions, we obtain the updated 
phase diagram shown in Fig. (22] The [12] state is favored 
when: 



b < a + 



t- 



M + 2b 
The [13] is favored when: 

e 



- at 4 - (4| 7 | + 3/3)t 4 (5.14) 
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FIG. 22: (Color online) Fourth order phase diagram. The 
new [13] state has width ~ t 4 . 



The columnar state is favored when: 



b> a 



Ad+2b 



at 4 + f3t 4 



(5.16) 



(5.15) 



On the [12]-[13] boundary, there is a degeneracy be- 
tween intermediate states where adjacent lines are sep- 
arated by either two or three columnar strips. On the 
[13]-columnar boundary, there is a degeneracy between 
states where adjacent staggered lines are separated by at 
least two columns. 



3. Higher orders and fine structure 

The picture of Fig. [22]will be further refined by consid- 
ering sixth order resonances and new phases will appear 
near the phase boundaries, in regions of width ~ t e , 
which is why they were missed at fourth order. The most 
immediate consequence will be the lifting of the degener- 
acy of the A3 states in favor of the state [13]. The latter 
state, in comparison with the other A3 states, is both 
stabilized maximally by the sixth order analog of Fig. [5TJ 
and, if we choose p > q, destabilized minimally by the 
sixth order analog of Fig. [20] The [13]-columnar phase 
boundary will open to reveal the [14] phase (Fig.[Sbl)^, in 
which the number of favorable long resonances, the sixth 
order analogs of Fig. [21k . is maximized and there are no 
repulsive contributions ( i. e. the sixth order analogs of 
Fig. [50] will be disconnected terms if the lines are more 
than three columnar strips apart). 

The argument may be applied iteratively at higher or- 
ders in perturbation theory. At 2n-th order, we may ask 
whether the [ln]-columnar boundary will open to reveal 
a new phase. The transition to less tilted states will be 
again be driven by processes that connect adjacent lines 
(Fig. |23|) . In the competitive states, adjacent lines are 
separated by at least n columnar strips so 2n-th order 
resonances connecting the lines must be "straight" . This 
means that the complicated high order processes, includ- 
ing "snake-like" fluctuations that break the staggered 
lines, will simply change the self-energy of a staggered 
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line and do not have any effect on the transition. The 
[In] phase will be destabilized by the process in Fig. [23b 
which will overwhelm the stabilizing effect of Fig.l2"3"a due 
to combinatorics. However, these repulsive processes will 
not contribute when the lines are separated by more than 
n columnar strips. Therefore, the [l,n+l] phase, which 
maximizes the number of the long 2n-th order resonances 
(Fig. I2"5b). will be stabilized in a region of width ~ t 2n 
between the [In] and columnar phases. Therefore, we 
obtain the phase diagram of Fig. [T5] 



(a) 



(b) 



(c) 



FIG. 23: (Color online) These are the 12-th order resonances 
which drive the transition between the [16] and [17] states. 
Process (a) selects the [16] state from the others in the A§ 
manifold, which were stabilized equally at 10th order. Process 
(b) destabilizes the [16] state near its columnar boundary in 
favor of the [17] state, which maximizes the number of most 
favorable resonances (c). 

This shows that states with arbitrarily long periods 
are stabilized without long range interactions or fine tun- 
ing (other than the requirements of perturbation theory) . 
The situation will be similar if we were to use the two- 
dimer Rokhsar-Kivclson resonance. In this case, there 
will be contributions from resonances occurring only in 
the columnar regions, which were "inert" in our calcula- 
tion. These processes will amount to self-energy correc- 
tions that just renormalize columnar energy scale b. Also, 
additional local terms (i.e. other than p and q) may be 
required to realize the very high-order states because ad- 
jacent lines will be able to interact via intermediate states 
other than the ones shown in Fig. [20] While we have not 
worked out the exact details of this case, we note that 
there are a finite number of such intermediates states so 
only a finite number of local terms will be required. In 
particular, we will not have to add longer terms at each 
order in perturbation theory. 

So far, we have concentrated on the boundary with the 
columnar phase but we may also ask whether a similar 



lifting may occur on the other phase boundaries. We 
consider the [11]-[12] boundary. Both of these phases are 
stabilized at second order and occupy regions of width ~ 
t 2 in the phase diagram (assuming \c—a\ 3> t 2 ). On their 
boundary, all states where staggered lines are separated 
by either one or two columnar strips are degenerate to 
second order. We investigate the effect of fourth order 
resonances on this boundary. 

We need to consider not only the resonances presented 
in section fVB 2l but also new fourth order processes which 
become available once we consider staggered lines that 
are one column apart. These are shown in Figs. [24] 
and [25] The resonances in section IV B 21 will stabilize 
(or destabilize - the sign is not important) each bound- 
ary state by an amount proportional to the number of 
its "[12] regions", i. e. columnar regions that are two 
columns wide and have staggered lines on their bound- 
ary, while the resonances in Fig. [24] will contribute an 
amount proportional to the number of "[11] regions", i. 
e. columnar regions that are one column wide. If these 
were the only available processes, the [11]-[12] bound- 
ary would be shifted by ~ t A but the degeneracy on the 
boundary would remain. 

The possibility of a new phase is determined by the 
resonances in Fig. [55] Both of these resonances have an 
overall stabilizing effect (the contributions to the energy 
are negative) and depend on whether a [11] region is ad- 
jacent to another [11] region (Fig. 125k) or a [12] region 
(Fig. [23b). Because c > a, resonance (b) is stronger 
than (a), since its intermediate state has lower energy, 
but requires a lower density of staggered lines. If the 
net contribution of resonance (a) wins, then the degen- 
eracy would be lifted in favor of the [11] state and the 
[11]- [12] boundary would be shifted again, but the de- 
generacy will only be between the two states. However, 
if c is made sufficiently large^i, the contribution of reso- 
nance (a) tends to zero while (b) approaches a constant 
because the intermediate state can occur without involv- 
ing c bonds (i. e. the fourth order process where the left 
cluster is flipped first and last). Therefore, there will be a 
new phase where resonance (b) is maximized. This state 
is the [11-12] state where the label refers to the repeating 
unit "one staggered strip followed by one columnar strip 
followed by one staggered strip followed by two columnar 
strips" . 

Continuing the line of thought, we may ask whether the 
[11-12]-[12] boundary will open at higher orders. Sixth 
order resonances will shift the boundary but there are 
no processes which break the degeneracy. However, at 
eighth order, there is a resonance which will favor a [11- 
(12) 2 ] phase (Fig. [26). 

We can, in principle, investigate whether the [ll-(12) n ] 
phases continue to appear when n is large and also 
whether the new boundaries themselves open to reveal 
even finer details. The same arguments will hold for 
all of the other boundaries in the [In] sequence. While 
the structure of our Hamiltonian allowed us to be def- 
inite regarding the [In] sequence, it is more difficult to 
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(a) 



2 flips 




(a) 



initial state 




FIG. 24: (Color online) These are fourth order processes avail- 
able between two staggered lines separated by a single colum- 
nar strip. 




(b) 



FIG. 26: (Color online) These are eighth order processes 
which can lead to a new phase between the [11-12] and [12] 
phases. Resonance (a) stabilizes the [11-12] phase while (b) 
stabilizes the [11-(12) 2 ] phase. In the limit where c is large, 
resonance (b) is preferable. The easiest way to see this is by 
setting c — oo. Then the energy term of resonance (a) and all 
of the wavefunction terms in resonances (a) and (b) will give 
zero because the intermediate states can not form without cre- 
ating c bonds. However, the intermediate state in resonance 
(b) does not involve c bonds and can be obtained without cre- 
ating c bonds (i. e. in an eighth order process where the first 
and last two actions involve flipping the cluster on the right 
and the first cluster in the middle). Therefore, the energy 
term of resonance (b) will give a stabilizing contribution. 



FIG. 25: (Color online) These are fourth order processes 
which can occur if a [11] region is next to (a) another [11] re- 
gion or (b) a [12] region. Both of these resonances contribute 
to the energy with a negative sign because the intermediate 
state involves two fewer repulsive c bonds than if the flipped 
clusters were farther apart. 



draw conclusions about the fine structure at high orders 
in perturbation theory because increasingly complicated 
resonances need to be accounted for. Most of these reso- 
nances will stabilize the unit cells of the states on either 
side of the boundary so the net effect will be to move the 
boundary. The more important terms, with respect to 
whether boundaries will open, are resonances associated 
with interfaces between regions with one or another unit 
cell. Even these terms can become complicated at very 
high orders in perturbation theory. However, our argu- 
ments suggest that arbitrarily complicated phases can, in 
principle, be stabilized by going to an appropriate range 
in parameter space and/or adding additional local in- 
teractions. Therefore, the most generic situation is an 
incomplete devil's staircase, as sketched in Fig. [TH 



C. Order of the transitions between the [In] phases 

The boundaries between different modulated states are 
generically first-order. Intuitively, this is not surprising 
because of the topological property that "protects" the 
states, namely that even with an infinite number of local 
dimer moves, it is impossible to go from one state to the 
other since the states arc in different topological sectors. 
Therefore, we would not expect a transition to be driven 
by a growing correlation length. 

Formally, the way we use to determine the order of the 
transitions that emerge in the system, is by calculating 
the first derivative of the ground-state energy on cither 
side of the phase boundaries that we founds For ex- 
ample, wc will treat explicitly the case of the transition 
between [12] — ► [13], even though the same line of argu- 
ment applies to the other boundaries as well. The ener- 
gies of the two states near the phase boundary, are given 
by Eg. I5.11[ Eq. l5.12l and the phase boundary is given by 
the following condition (as can be seen in Eg. I5.14|) : 

& = a + ^-^-ai 4 -(4|7|+3/3)* 4 (5.17) 

Let's consider the case where we approach the boundary 
from the [12] side varying the variable t but keeping a, b 
constant. In the phase diagram Fig. [T31 we 'move' ver- 
tically down. The reason for choosing this direction is 
just clarity. Let's call the point of the phase boundary 
where our path crosses, A, and the critical value of t, t c 
(coming from the solution of Eq. 15.171 for fixed values of 
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a, 6, d). 

The energies of the two states at the phase boundary 
are exactly equal. Their derivatives are: 



dE [12] 



dt 



= -( 
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By Eqs. I5T71 I5J81 and [5T9l we have 
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<0 to order 0(t 2 ) 



The derivatives are not equal along the phase bound- 
ary so the transition is discontinuous (first-order). It is 
clear that all the phase transitions we found will also be 
discontinuous because the above discontinuity comes ex- 
actly from the contributions leading to the phase bound- 
ary's presence. 



VI. CONNECTIONS WITH FRUSTRATED 
ISING MODELS 

A natural question to ask is whether the staircase pre- 
sented above has any connection to the staircase of the 
3D ANNNI model or the quantum analogs discussed in 
Ref. [5l|. One of the main differences of the present work 
is the non-perturbative inclusion of frustration by con- 
sidering hard-core dimers as the fundamental degrees of 
freedom. In this sense, the present staircase differs from 
previous work similarly to how the fully frustrated Ising 
model differs from the conventional Ising model. It is in- 
structive to consider the mapping between dimer cover- 
ings and configurations of the fully frustrated Ising model 
on the square lattice (FFSI) in more detail. 

The FFSI model can be described in terms of Ising de- 
grees of freedom living on the square lattice. The main 
difference with the usual ferromagnetic Ising model is the 
following: Even though in the x-direction. all the bonds 
are ferromagnetic, in the y-direction there are alternating 
ferromagnetic (antiferromagnetic) lines where ferromag- 
netic (antiferromagnetic) vertical bonds live (we consider 
that the absolute values of the couplings of all bonds are 
equal). In this way, the product of bonds on a single 
plaquette is always —1 (three ferromagnetic bonds per 
plaquette) and therefore the ground-state of the system 
cannot be just the ferromagnetic one. In fact, by map- 
ping each " unsatisfied" bond to a dimer living on the dual 
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FIG. 27: (Color online) One of the four possible [11] config- 
uration in terms of Ising spins on the fully frustrated square 
lattice. The hard-core dimer constraint corresponds to the re- 
quirement that the FFSI ground state has one "unsatisfied" 
bond per plaquette. The columnar-dimer regions correspond 
to ferromagnetic-Ising domains. The staggered-dimer strips 
correspond to the Ising domain walls, depicted by the red- 
colored (dashed) curved lines, which separate ferromagnetic 
domains of different orientation. Clearly, in the other equiva- 
lent configurations, even though they share the same principle 
of domain- wall competition, the separated regions are not fer- 
romagnetic but one of several metamagnetic choices^. 



sublattice, we find that each of the degenerate ground- 
state configuration maps to a hard-core dimer configura- 
tion on the square lattice (see Fig|27|^. 

A typical [In] configuration in the dimer language we 
used, as clearly depicted in Fig. [27] for one of the four 
equivalent dimer structures, under 7r/2 rotations and sub- 
lattice shifts, can be seen in the FFSI picture as ferro- 
magnetic stripes of length An in the one direction and 
infinite in the other, separated by antiferromagnetic do- 
main walls. In this way, these ordered states clearly re- 
semble the modulated phases of the ANNNI model in 
two dimensions. The other three equivalent dimer struc- 
tures map again to periodic domains in the Ising model, 
but with more complicated metamagnetic structure. The 
reason for this seemingly large complexity has to do with 
the fact that the possible equilibrium configurations have 
to satisfy the FFSI constraint. 

As far as the interactions are concerned, we have the 
following correspondences: The three-dimer resonance 
term we used in our construction corresponds to a two 
neighboring-spins flip process which should, however, re- 
spect the FFSI constraint. We should note that the usual 
single-plaquette resonance move maps to the single-spin 
flip which is the same as the Ising transverse field usually 
considered. The a and 6-tcrms correspond to domain- 
wall energies in the FFSI. Both of them have an one- 
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to-one mapping to three-spin interactions but these in- 
teractions are also anisotropic (they depend on the dis- 
tribution of the Ising bonds which we described). The 
additional interactions that we added to the system, so 
that to extensively study it, correspond clearly to com- 
plicated multi-spin interactions. 



VII. DISCUSSION 

There are reasons to be optimistic that these ideas ap- 
ply more generally. For example, as mentioned earlier, we 
expect that with suitable modification of Eq. 15.11 a simi- 
lar phase diagram may be obtained for a wide variety of 
off-diagonal resonance terms, including the familiar two- 
dimer resonance of Rokhsar-Kivclson. This is because 
the perturbation theory is structured so that at any or- 
der, most of the nontrivial resonances amount to self- 
energy corrections and the resonances driving the tran- 
sitions are comparatively simpler. The three-dimer reso- 
nance of Eq. l4.2l is analytically convenient as its action is 
confined to the domain wall boundaries. The two-dimcr 
resonance would involve more complicated bookkeeping 
since we also need to account for internal fluctuations of 
the columnar regions. 

Another reason to expect these ideas to hold more gen- 
erally is the qualitative similarity of this approach to the 
field theoretic arguments in Refs. and |29|. In those 
studies, the following action, in the notation of Ref. l29l . 
was used to describe the tilting transition in the Rokhsar- 
Kivelson quantum dimer model on the honeycomb lattice 
(the square lattice is similar but with some added sub- 
tleties - see Ref. [H): 

S = l -{d T hf + \p2(Vhf + l - Pi ^ 2 hf 

+ 93(d x h)(^d x h - ^-d y h){^d x h + ^-d y h) 
+ g A \Vh-Vh\ 2 + ... (7.1) 

where ". . ." includes terms that are irrelevant to the 
present discussion (though maybe not strictly "irrele- 
vant" in the RG sense). In this expression, h is a coarse 
grained version of the height field (Fig. [5]) and the first 
line of Eq. 17.11 describes the tilting transition at the RK 
points, which corresponds to p2 — gz = g± = 0. If 
<?3 < 0, the system favors tilted states and is similar to 
our parameter a — b. However, g<± prevents the tilt from 
taking its maximal value and in this sense, is similar to 
our terms c and d. The existence of the devil's stair- 
case in Ref. Hi was established by tuning g 3 and 34 so 
as to stabilize an intermediate tilt and then to study the 
fluctuations about this state. The staircase arose from 
a competition between these quantum fluctuations, anal- 
ogous to our term t, and lattice interactions, (roughly) 
analogous to our terms c, d, p, and q. 

Another sense in which our calculation is similar to 
Ref. Hi may be seen by heuristically considering the effect 



of doping the model. In particular, consider replacing one 
of the dimers in a staggered strip with two monomers. If 
we then separate the monomers in the direction parallel 
to the stripe, a string of columnar bonds will be created. 
If the staggered and columnar bonds were degenerate, 
then this would cost no energy in addition to the cost of 
creating the monomers in the first place so the monomers 
would be deconfined. However, in the [In] phase, the 
staggered bonds are slightly favored so the energy cost 
E of separating the monomers by a distance R would be 
E ~ Rt 2n . Hence, the commensurate phases seen in our 
model are confining with a confinement length that be- 
comes arbitrarily large for the high-order structures that 
appear very close to the columnar phase boundary. This 
is qualitatively similar to the "Cantor dcconfincmcnt" 
scenario proposed in Ref. I29I . 

However, there are ways in which our calculation is 
qualitatively different from the above. Our calculation 
takes place in the limit of "strong-coupling" where t 
is small compared with other terms but influences the 
phase diagram nonetheless because the stronger terms 
are competing. In contrast, the RK point of a quantum 
dimer model, by definition, occurs in a regime of param- 
eter space where quantum fluctuations are comparable in 
strength to the interactions. The field theoretic predic- 
tion requires g% and <?4 being nonzero so does not liter- 
ally apply at the RK point cither but, by self-consistency, 
should apply somewhat "near" it. We may speculate that 
the tilted states being predicted by the field theory arc 
large t continuations of states that emerge in the strong 
coupling limit far from the RK point. However, we reem- 
phasize that our calculation is reliable only in the limit 
of small t and we can not be certain which (if any) of our 
striped phases survive at larger t. Another issue is that 
the phase diagram near the RK point depends strongly 
on the lattice geometry and the prediction of a devil's 
staircase in Ref. Hi is for bipartite lattices. In contrast, 
lattice symmetry does not play an obvious role in the 
present work and it is likely that these ideas can be made 
to apply on more general lattices. 

It is also likely that this calculation can be made to 
work in the strong-coupling limits of other frustrated 
models, for example vertex models^! and even in higher 
dimensions. For example, mappings similar to those dis- 
cussed in Ref. [H may be used to construct an SU(2)- 
invariant spin model on a decorated lattice that displays 
the same phases. A more interesting direction would 
be to study the strong-coupling limits of more physi- 
cal models, for example the Emery model of high T c 
superconductivity^ which also shows an affinity for ne- 
matic ground states^. It would also be interesting to see 
whether nematicity is essential, i. e. whether other types 
of phase separation can occur in a purely local model 
through effective long range interactions that arise, as 
in the present calculation, from the interplay of kinetic 
energy and frustration. 
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APPENDIX A: DETAILS OF THE SECOND 
ORDER CALCULATION 

In this section, we explicitly work out the details of the 
second order calculation, including the effect of having 
c > a. We show that the second order phase diagram is 
qualitatively similar to the c = a case presented in the 
main text, provided that c is not too large (the precise 
condition is obtained below). Therefore, fine-tuning to 
c = a is purely a matter of convenience. 

With reference to Fig. [THl we may calculate the differ- 
ence in energy between the excited and initial states for 
each case. The unperturbed energies e\ n , i G {1, 2, 3}, of 
the initial states in Fig. HSU. fTBT-2. and [TBI- 3 are given 
by Eq. 15.21 and are degenerate when a = 6. To calculate 
the unperturbed energies, ef* of the excited states, we 
need to examine the interaction energies present in the 
excited states which are absent in the initial states and 
vice versa. These are shown in Fig. by the red and 
blue arrows. Using the figure, we obtain: 



ex 
e l 



(Ad -2a) - (-2a- 26) 

4d + 26 (Al) 



= (Ad + 2c - 2a - 2a) - (-2a - 2b) 

= Ad + 2b + 2c - 2a (A2) 

ef-ef = (6d + p - 2a - 6) - (-2a - 26) 

= 6d + p + b (A3) 

These are the three possible energy denominators which 
may enter Eg. 15. 41 We may classify a staggered line based 
on the dimcr arrangement to its immediate right, cor- 
responding to cases (1), (2), and (3) in Fig. [L6j Each 

.2 

staggered line contributes — N y e J_ in . The factor N y 



is because there are ^ possible rightward resonances of 
each line and also ^ leftward resonances of the line to 



its right, which enter with the same weight. Putting ev- 
erything together, we can update Eq. 15.21 



E(N sU N s2 ,N s3 ,N c ) = ~b^^ + (b-a)N y N s 



- N y N sl 



NyN s2 



NyN s3 



4d+26 



i 1 



Ad + 26 + 2c - 2a 
t 2 



6d + p + b 
+ 0(t 4 ) (A4) 

where N S i is the number of staggered lines of type (i) and 
N s = N s i + N s2 + N S 3 is the total number of staggered 
lines. While the zcrocth order term depends only on the 
number of staggered lines, the 0(t 2 ) piece depends on 
their distribution and relative orientations. 

Eq. lA4l shows that type (1) staggered lines will be sta- 
bilized more by the perturbation than type (2) lines. 
However, the zeroeth order term depends only on the 
total number of lines and type (2) lines permit a denser 
packing of lines. We temporarily ignore type (3) lines. 
Depending on 6 — a, the system will favor either the [11] 
state, the columnar state, or the states which maximize 
the number of type (1) lines, which (at second order) 
are analogous to the maximally staggered configurations 
in Fig. [7] except staggered lines are now separated by 
two columns. We denote the latter collection of states 
with the label A 2 (i. e. alternating states where staggered 
strips alternate with two columns of columnar dimcrs). 
To determine these boundaries, we first write down the 
energies of these three states to second order in t using 
Eq.EU 

N y N x N„N r . 



E[n] — —b- 



Ea 2 



(b-a 



f- 



Ad + 2b + 2(c-a)' ) 
(A5) 



NyN X 



NyN x 

— (6 - a 



NyN x 



Ad + 2b' 



(A6) 



(A7) 



where we have used N s = ^ for the [11] state, N s = ^ L 
for the A 2 states, and N s = for the columnar state. 
Comparing these expressions, we obtain the following 
boundaries. The [11] state is favored when: 



6 < a 



3i 2 4d + 26 



Ad + 2b 2(4d+26) v 4d + 26 + 2(c-a) 

(A8) 

The A 2 states are favored when: 

t 2 it 2 Ad + 2b 

+ 4d+26 ~ 2(4d+26) ( ~ 4d + 26 + 2(c - a) ' < 



< a + 



Ad + 2b 
(A9) 
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and the columnar state is favored when: 



b> 



4d+2b 



(A10) 



If we require that (c — a) is small compared to 2d + b 
(recall that d > a by assumption), then we see that the 
region where the A2 state is preferred is a small region 
within the phase boundary of Fig. [TUJ The width of this 
region tends to zero as c — ► a. A cartoon of this case 
is given in Fig. 1171 On the ^-columnar boundary, the 
A 2 states are degenerate with the columnar state and 
any intermediate state where consecutive staggered lines 
arc separated by at least two columns. Similarly, on the 
[ll.J-.A2 boundary, the [11] states are degenerate with the 
A2 states and any intermediate states where staggered 
lines of different (same) orientation are separated by two 
(either one or two) columns. 

The collection of A2 states are degenerate to order t 2 . 
As the degeneracy of the A2 states will be partially lifted 
at fourth order in the perturbation theory, the c > a 
case adds complexity without changing the phase dia- 
gram qualitatively. Therefore, purely for convenience, 
we assume c = a throughout the main text. 

As c is made larger, the width of the A2 region in- 
creases but when c > d + a + ^j^, then according to 
Eq. IA41 we need to consider type (3) lines. In this case, 
the tilted state will no longer be favored and, in fact, all 
of the steps of the staircase will lie in the zero winding 
number sector. 



APPENDIX B: CANCELLATION OF 
DISCONNECTED RESONANCES 

We now demonstrate the cancellation of disconnected 
terms that appear at fourth order in the perturbation 
theory (Eq. I5.8|) . Because our Hamiltonian is local, there 
are linked cluster theorems which ensure that this can- 
cellation occurs at any order in the perturbation theory 
so the contribution to the energy is always extensive as 
it should be. 

Regarding the situation of Fig.ll8| we may write down 
all fourth order terms involving the excited states which 
we have labelled 1 and 2. The energy numerator will 
contribute a number of terms, each having energy: 



E 



-t 4 



VnmVmlVlkVkn 



(e m - e n )(e; - e n )(efc - e n ) 

-t 4 



(4d + 2&) 2 (2(4d+26)) 



(Bl) 



which follows because the energy of the excited states 
(relative to the initial state) are ei — eo = (-2 — eo = 4d+2b 
and ei2 - eo = 2(4d + 26). That ei 2 — eo = 2(ei - e ) is 
precisely because clusters 1 and 2 are disconnected. This 
would not be the case if we had a long range interaction. 
From Fig. [TH we readily see that there will be four such 
terms since there are four ways to connect the initial 
state (0) to the excited state (12) and then back to itself. 
Therefore, the full potentially pathological contribution 
of the energy numerator is: 



E n 



t: 1 



t 4 



"(4d + 2&) 2 (2(4d + 26)) 



(4d- 



2b) 3 
(B2) 

It is evident that there are of order N 2 of these terms 
because the choice of clusters 1 and 2 were arbitrary. 
The energy denominator (wavefunction normalization) 
will contribute terms, each having energy: 



E = f 



VnrnVrnnVnlVln £ 

(e m -e„) 2 (e,-e„) = (4d+26)3 



(B3) 



and only resonances from the initial state to excited 
states 1 and 2 are involved. Because the case where 
to = I gives an extensive contribution to the energy, we 
will be concerned with the case where to and I in the 
above equation are different. There are two such terms 
because either to or I can correspond to excited state 1 
and the other to state 2. Therefore, the full potentially 
pathological contribution of the energy denominator is: 



E, 



den 



{Ad + 2bf 



(B4) 



We see explicitly that the non-extensive contributions of 
E num and Eden precisely cancel. Clearly this will be case 
for any choice of disconnected clusters 1 and 2 so we have 
shown that, to fourth order in the perturbation theory, 
the energy correction is extensive, as it physically should 
be. 
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